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Pairing between fermions that attract each 
other, reveal itself to the macroscopic world in 
the form of superfluidity. Since the discovery 
of fermionic superfluidity [l|, intense search has 
been going on to find various unconventional 
forms of fermion pairing (2| as well as to in- 
crease the transition temperature fS\. Here, we 
show that a two dimensional mixture of single- 
component fermions and dipolar bosons allows 
to reach experimentally feasible superfiuid transi- 
tion temperatures for non-standard pairing sym- 
metries. Excitations in these superfluids are any- 
onic and their statistics depends on the order of 
their permutations, i.e is non-Abelian. Our re- 
sults provide for the first time an example of a 
highly tunable system which exhibits various kind 
of pairing symmetry and high transition temper- 
ature. Additionally, they provide a playground 
to observe anyonic excitations and their braiding 
properties. 

Attraction between fermionic particles favours pairing 
of fermions resulting in superfluidity of the system. The 
paired fermions, known as Cooper pairs, can have differ- 
ent kind of internal symmetries. The common ones found 
in nature have s-wave and c?-wave like internal structure 
and conserves parity and time reversal symmetry. Also, 
Cooper pairs with chiral Px + ipy-wave internal struc- 
ture have been proposed for the observed superfluidity 
of electrons in Strontium Ruthnate [3]. This kind of 
pairing breaks the time-reversal symmetry. Spin-less chi- 
ral p-wave superfluid state has formal resemblance with 
the "Pfaffian" state proposed in relation to the fractional 
quantum Hall state with filling factor 5/2 When 
confined in a two dimensional geometry, excitations in 
the chiral p-wave superfiuid become non-abelian anyons. 
Anyons are particles living in a two dimensional plane 
that under exchange do not behave as bosons or fermions. 
In some cases, exchange of two su ch p articles depends on 
the order of the exchange [1, H, [ij . Particles obeying 
such laws are called non-Abelian anyons. Apart from 
fundamental interest on the occurrence of such particles, 
non-Abelian anyons find remarkable applications in the 
field of quantum information for quantum memories and 
fault-tolerant quantum computation flU . Recently it has 
been shown that quasi-particles living in stable vortex ex- 
citations of chiral two-dimensional p-wave spinless super- 
fiuids obey non-Abelian statistics jH, [l^ . Using p-wave 
Feshbach resonances in fermionic ultracold atoms, such 
superfiuids can be realized in principle, but this proce- 



dure is very difficult because of non-clastic loss processes 

Bose-Fcrmi mixtures arc another candidate for creat- 
ing superfluidity in fermions via boson mediated inter- 
actions and have formal resemblance with phonon medi- 
ated superconductivity in metals It was found, how- 
ever, when the bosons and single-component fermions 
are completely mixed, the maximum possible transition 
temperature is of the order of 10~^Tp for p-wave pair- 
ing, where Tj? is the Fermi Temperature [l5|. Any at- 
tempt to increase the transition temperature by increas- 
ing the boson-fermion interaction strength or fermionic 
density results in phase separation between the mixture. 
Hence it is experimentally hardly possible to study this 
phenomenon for unconventional superfluidity in ultracold 
atomic systems. 

Here, we show a possible way of overcoming these dif- 
ficulties. We study the property of superfiuidity in Bose- 
Fermi mixtures, where bosons are interacting via long- 
ranged dipolar interactions. We show that the transition 
temperature for p-wave superfluidity can become compa- 
rable to the Fermi energy. We find that other more exotic 
Cooper pairs with /- and /i-wave internal symmetries are 
possible in certain range of Fermi energies without bosons 
and fermions separating. To the best of our knowledge, 
for the first time a system is proposed where conventional 
pairing mechanism gives rise to different exotic internal 
structures of the Cooper pairs with strong interaction in 
respective angular momentum channels. In addition, we 
study the excitations in chiral states of the odd-wave su- 
perfiuids and point out their non-Abelian anyonic nature. 

Experimentally, an available bosonic species, where 
prominent dipolar interaction can be achieved using Fes- 
hbach resonance is Cr^^ [13, [H] . Another route towards 
achieving dipolar condensate is to experimentally realize 
quantum degenerate heteronuclear molecules [19| which 
have permanent electric moment. Thus in the near fu- 
ture a quantum degenerate mixture of dipolar bosons and 
fermions will be achievable experimentally. 



I. SYSTEM 

A. Dipolar Bose Condensate 

Our system, as sketched in Fig[T^, consists of dipolar 
bosons mixed with single component fermions confined 
in a quasi-two dimensional geometry by a harmonic po- 
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tential with frequency tOz and oscillator length Iq. Here 
we present a brief overview of dipolar condensates with 
a focus on our present problem (for details sec meth- 
ods). First, we assume that the bosons are polarized 
along the z direction. The dipolar interaction reads 
Vdd = 47rgdd ^^j^2 ^ _ jj^ momentum space, where 
gdd is the dipole-dipolc interaction strength. For atoms 

Odd = MoA^Ti 

where fj,m and are the magnetic moment of the atoms 
and the electric dipole moment of the molecules, respec- 
tively. We assume that the z dependance of bosonic 
wavefunction is given by a Thomas-Fermi profile with 
radius Rz- By integrating out the z direction, the effec- 



fcrmion interaction becomes function of momentum. In- 
tegrating out the bosonic degree of freedom results in 
effective interaction between the fermions llal. 



Vph((7l,t^) = 



95, 
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Ubq'^/mb 
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where fc^ 



/47r, and for dipolar molecules gdd = Me/47i-eo ^^'^^^ ^ .^^'K^ " '^)' momentum exchange 

between the interacting particles along the Fermi surface, 
and a is a number which depends on Rz/^f (see Meth- 
ods). This form of interaction has formal resemblance 
to superconductivity in metals. Due to the momentum 
dependance, if attractive, Vpii{qj_,uj) makes the Fermi 
surface unstable against formation of Cooper pairs with 
higher internal symmetries. Assuming momentum trans- 
fer around Fermi momentum kj, in two-dimensions we 
can expand Eq. ^ as 



tivc dipolar interaction takes the form V^ff 



V(/cj 



kl + ky and V{k±) denotes the shape of the 
interaction. Next, we define a dimensionless dipolar in- 
teraction strength g^d = 8Trmbgddnbio/5h'^ which will be 
used later, where £0 is the ground state oscillator length. 
V{k±) is repulsive for small momentum and attractive 
in the high momentum limit. At low temperature, the 
ground state of the bosons is a condensate with fluctua- 
tions around this state. Generally, the spectrum of these 
excitations, denoted by fl(kj_), can be divided into two 
parts: i) ~ fcj_, phonon spectrum for small momenta and 
ii) ~ /c^, free-particle like spectrum for higher momenta 
[20j. For the dipolar condensates and g^d greater than 
a critical value, Q{kj_) has a minimum at momentum 
^0 where fcp is in intermediate momentum regime 
as shown in Fig. 2b. Following Landau, the excitations 
around the minimum are called "rotons" . With increas- 
ing g^d the excitation energy at fco decreases and even- 
tually vanishes for a critical particle density as shown in 
Fig. 2b. When the particle density exceeds that critical 
value, the excitation energy becomes imaginary at finite 
momentum and the condensate becomes unstable. The 
"roton" part of the spectrum is absent in condensates 
with contact interaction. 



B. 



Effective interaction 



In this section we discuss the effect of bosons on 
fermions. The fermions are interacting with the bosons 
via short ranged contact interaction of strength g^f. 
The fermionic density dependance on z is given by a 
normalized Gaussian with width if. Wc also define 
the dimensionality parameter for the fermions as 77 = 
fi/hiOz, where ^, is the chemical potential for the fermions. 
For < 1 the fermionic system is quasi-two dimensional, 
and we can assume £f ~ y^H/mfOJz, i.e the width is 
given by the oscillator length of the trapping potential. 
77 > 1 corresponds to three dimensional fermionic sys- 
tem and ^/ is given by minimizing the Thomas- Fermi en- 
ergy functional, if ^ Tf'^'^{h/mfUJzY^'^ ■ The fluctuations 
in fermionic density couple to the density fluctuations 
present in the Bose condensate. Due to the momentum 
dependance of the bosonic excitations, the condensate- 
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where the dimensionless effective interaction between the 
fermions in angular momentum channel m is given by 
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and A^o = mf/ir'^h^ is the two-dimensional density of 
states for the fermions. In the present paper. Am > de- 
notes attractive interaction. As we are considering single 
component fermions, pairing will occur in the odd an- 
gular momentum channels. The effective dipole-dipole 
interaction between the bosons is given by the density 
dependent term g^d- The expression for effective inter- 
action obviously makes sense as long as the excitation 
frequencies of the bosons are real. To look into the prop- 
erties of effective interaction, we consider first the in- 
teraction strength in p-wave channel Ai. As shown in 
Fig. 2b, increasing g^d results in smaller roton gap, and 
when (73d ~ 3.61, the roton minimum touches the zero 
energy axis. Subsequently the effective fermionic inter- 
action Ai increases and diverges as (73d — > 3.61, as shown 
in Fig. 2a. Depending on the dimensionality 77, the rate 
of divergence changes and the interaction can become re- 
pulsive. This Feshbach resonances-like characteristics are 
the novel feature of this system. 

Next wc look into the variation of Am as a function of 
77 as depicted in Fig. [3] Here for concrctencss we con- 
sider the particular case corresponding to a Chromium- 
Potassium mixture with nib = 52a. m.u. and nif = 40 
a.m.u.. Using these masses, the interaction strengths in 
the channels m = l(p-wave),= 3(f-wave),= 5(h-wave) 
have been plotted in Fig. 3. As 77 — > 1, interaction in 
the p-wave channel becomes predominant. This trend 
is also followed when we are inside the three dimensional 
limit with ?/ > 1. But, with decreasing dimensionality we 
find that the predominant interaction channel changes 
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(a) (b) 

FIG. 1; a) Schematic image of a two-dimensional mixture of dipolar bosons (semi-transparent spheres with arrows) and fermions 
(blue spheres), b) The interaction between the bosons and fermions induces interaction between the fermions. This results in 
superfluidity of the fermions with various angular momenta as inducated by the arrows around the fermions. 




(a) (b) 

FIG. 2: a) Effective interaction strength in the p-wave channel Ai as a function of boson-boson effective interaction strength 
<73d. The green, blue and black lines correspond to 77 = 0.50, .60, .90 respectively, b) The excitation spectrum i^{k±) of the two 
dimensional Bose gas plotted as a function of k±. The blue, black and red lines correspond to gsd = 2.0, 3.0, 3.61, respectively. 



surprisingly from m = 1 to m = 5, the /i-wave chan- 
nel. Then A5 goes through a maximum attraction around 
77 ~ .6 and the interaction strength in the p-wave channel 
becomes repulsive. In this region, the repulsive p-wave 
interaction will renormalize the mass of the fermions and 
the supcrfluid instability will be due to m = 5 chan- 
nel. By decreasing 77 further, both Ai,5 become repulsive 
whereas the /-wave channel with m = 3 becomes attrac- 
tive. Decreasing rj further results in a situation, where all 
the channel have negligible interaction. The interaction 
strength at individual channel can be increase further by 
increasing the bosonic density closer to the critical 
value. 

Another obvious contribution to the strength of the 
interaction is the boson-fermion contact interaction 
strength gi,f. The condition for mechanical stability of 
the Bosc-Fcrmi mixtures reads [la 



-a^ < 1, 



8^.9ddV(0)4 

which constraints the magnitude of g\,f. From Fig. 2a 



and Fig. 3 it is clear that the maximum attractive inter- 
actions in each channel can be Ai,A2,A3 > 1 which are 
within the strong-coupling regime. Consequently maxi- 
mum transition temperatures possible in all three chan- 
nels are close to the fermi temperature. At this point 
let us compare this to the case of usual boson-fermion 
mixture, where highest interaction strength possible is 
Ai ~ .1 p^. [l6j. The maximum transition temperature 
possible is for p-wave symmetry and Tc is of the order 
~ exp(— ~ 10~^^. Note, then we have an im- 
provement in Tc of the order of 10^. 



II. CHIRALITY AND NON-ABELIAN ANYONS 



By considering the superconducting gap equation at 
low temperature, the gap is maximum when the order 
parameter breaks time reversal symmetry 0] . From now 
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on we assume that the order parameters are denoted by 



FIG. 3: Figure of effective interaction strength in the channel 
m = l(black hne),= 3(green hne),= 5(blue line) as a func- 
tion of the fermion dimensionality parameter rj. We fixed 
the bosonic effective interaction strength at pad = 3.1. The 
negative value of the interaction reflects it's repulsive nature. 



A,„ = Ao(r) 



k 



where = kcoaO, ky = fcsin0 and Ao(r) is the cen- 
ter of mass amplitude of the Cooper pairs with r be- 
ing the center of mass coordinate of the pair. For vor- 
tex state Ao(r) is defined as: i) Ao(r) = 0,r < ^ and 
ii) Ao(r) = Ag exp(i(/)), r > ^, where r = +2/^, 
tan0 = y/x. ^ is the size of the core of the vortex. 
The vortex state of the p-wave superfiuids always has a 
zero-energy bound quasi-particle state 0, [2^ . [23| . Now 
we discuss the asymptotic solutions for the zero-energy 
bound state for /- and ft,-wavc order parameters. The 
quasi-particle states in a single vortex can be found in 
the limit of large distance from the vortex core by solv- 
ing the Bogoliubov-DeGennes equation, 



d id 
dr r d(f) 

d id 

dr r d(j) 



(3) 



where E is the energy of the quasi-particles denoted by 
Um,Vm- We particularly look for zero energy solutions 



with bounded u„ 



with the property u„ 



The asymptotic solution of Eq. ([3]) with different orbital 
symmetries read, 



the vortices than the case for p-wave for similar values of 
gap and fermi energy. 
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The zero-energy solution for each odd-wave parameter 
corresponds to different angular momentum channel of 
the quasi-particles inside a vortex core. These results 
can also be carried out by applying the "index theo- 
rem" [l^l. For temperature smaller than the energy 
gap Ag//i, only the zero energy mode is occupied. The 
quasi-particle operator in that situation is written as 
7m = f (Pr(um(r)c^ (r) + Vm(r)c(r)), which acts as Majo- 
rana fermion fl^ . |22 | . 7„ obeys non-Abelian statistics 
and can be used for quantum computing [l^ . In order to 
perform quantum computational task, existence of sev- 
eral well separated vortices is necessary. From Eq. ^ we 
find that the condition for non-overlapped states can be 
achieved in m = 3 states with higher distance between 



IV. METHODS 

First we discuss the Hamiltonian describing the 
bosonic system which is homogeneous in the x — y plane 
and trapped in the z direction by a harmonic potential 
with frequency lOz- In the considered regime of parame- 
ters the bosonic density nb{x,y, z) is given by 



3n6 / z2 



where the Thomas-Fermi radius Rz is determined vari- 
ationally. By minimizing the mean field energy of the 
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BosG condensate within Thomas-Fermi regime we find 
that Rz/iQ = (Sgsd/S)^/"^. After integrating over z de- 
pendance of the density profile of bosons, the dipolar 
interaction takes the form Vcs = ^^^^{ki_) where 



1 

8 
15 



6(l + fci)cxp(-2fci) + 6 



where 



and 



are 



fcrmionic creation and destruction operator. ef{k±) = 

fc_L /2rnf is the dispersion energy of the fcrmions with 
mass nif. The density profile of fermions along the z 
direction is approximated by a gaussian with width i / . 

Including the fluctuations in the Bose condensate in 
the X — y plane, the condensate-fermion interaction 
Hamiltonian can be written as 



and k± = k±Rz and g is the contact interaction between 
the bosons which is assumed to be negligible in our case. 
Subsequently we write the Hamiltonian of the dipolar 
bosons in the condensed phase, Hj, = J2k]_ ^{k±)b''^-' 6^.- , 
where fe^- and 6^,- are Bugoliubov operators. The exci- 
tation spectrum is given in the units of trap frequency, 

n'ik^io) = + ff3d|^v [k^io?. 

Next, we consider the Hamiltonian describing the 
fermions and the boson-fermion interaction. Ki- 
netic energy for the single component non-interaction 
fermions, is characterized by the Hamiltonian Hf — 



k±,q1 



where the momentum dependent coupling constant is 



iven by j{k±) = y 2nbeh{k±) /Q,{k±) and 
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with crf(..) being the error function. 
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